The computational complexity of PEPS 
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We determine the computational power of preparing Projected Entangled Pair States (PEPS), as 
well as the complexity of classically simulating them, and generally the complexity of contracting 
tensor networks. While creating PEPS allows to solve PP problems, the latter two tasks are both 
proven to be #P-complete. We further show how PEPS can be used to approximate ground states 
of gapped Hamiltonians, and that creating them is easier than creating arbitrary PEPS. The main 
tool for our proofs is a duality between PEPS and postselection which allows to use existing results 
from quantum compexity. 



I. INTRODUCTION 

Computing the properties of correlated quantum 
many-body systems is a central task in many fields in 
physics. Its complexity stems mainly from the large 
dimension of the Hilbcrt space which grows exponen- 
tially in the system size. In the last decades, the Den- 
sity Matrix Renormalization Group (DMRG) method 
has proven extremely successful in the description of 
one-dimensional phenomena Recently, it has been 
shown that from the perspective of quantum information, 
DMRG can be described as a variational method over the 
class of Matrix Product States (MPS) 0. MPS struc- 
ture the state space into a hierarchy of states with poly- 
nomial description complexity [||, and it turns out that 
already the lowest levels of this hierarchy approximate 
many physical states of interest extremely well. MPS 
have a natural extension to two and higher dimensional 
lattices, called Projected Entangled Pair States (PEPS), 
which also have an efficient description and are promis- 
ing candidates for variational methods in higher dimen- 
sions [If. It has been shown that MPS can be created ef- 
ficiently by a quantum computer and that they also 
can be simulated efficiently classically @. In contrast, 
in two or more dimensions it seems to be hard to create 
arbitrary PEPS, as well as to classically compute expec- 
tation values. In fact, it has been shown that there exist 
2D PEPS which encode solutions to NP-completc prob- 
lems Q, thus posing lower bounds on their complexity 
and computational power. 

In the present work, we determine both the power of 
creating PEPS and the complexity of classically simu- 
lating them. We investigate which kind of problems we 
could solve if we had a way to efficiently create PEPS, and 
find that these are exactly the problems in the complexity 
class PP (deciding whether a boolean formula has more 
satisfying than non-satisfying assingments). Second, we 
show that classically computing local expectation values 
on PEPS is a #P-complete problem (counting the sat- 
isfying assignments of a boolean formula). This result 
can be extended to the contraction of arbitrary tensor 
networks, which turns out to be #P-comlete as well. 

The main tool in our proofs is a duality between PEPS 
and postselection, which permits to use existing results 



from quantum complexity [8|: any PEPS can be created 
by a postselected quantum circuit, and any output of 
such a circuit can be written as a PEPS. We also ap- 
ply this duality to show that ground states of gapped 
local Hamiltonians in D dimensions can be efficiently 
approximated by the boundary of a D + 1-dimcnsional 
PEPS. Finally, we compare the power of creating PEPS 
to the power of creating ground states of local Hamilto- 
nians. While in general they are equally hard, we find 
that when restricting to gapped Hamiltonians, creating 
ground states becomes easier: it is in the weaker class 
QMA, the quantum analogue of NP. 



II. PEPS AND POSTSELECTION 

We start by recalling the definition of PEPS [§]. Con- 
sider an arbitrary undirected graph where each of the 
vertices corresponds to a quantum system (a spin) of 
Hilbert space dimension d. A PEPS on these N spins 
is constructed by placing as many virtual spins of di- 
mension D on each vertex as there are adjacent edges. 
Along each edge, these virtual spins form maximally cn- 
tanged states X)t=i I*) I*)- The physical spins are now 
obtained from the virtual ones by applying a linear map 
PH : C D <g> . . . <g> C D -> C d at each vertex v. For the sake 
of readability, we will mostly supress the dependence of 
P on v. The graph underlying the PEPS will usually be 
chosen according to the physical setup, typically a two 
or higher dimensional lattice. 

Let us now turn to postselected quantum circuits 
Roughly speaking, postselection means that we can mea- 
sure a qubit with the promise of obtaining a certain 
outcome. More precisely, the postselected circuits we 
consider start from the 1 - - - 0) state, perform a se- 
quence of unitary one- and two-qubit gates, and post- 
select on the first qubit being |0). Thereby, the state 
a|O)|0o)+/?|l}|</ ) i} is projected onto the state |0o), which 
is the state created by the postselected quantum cir- 
cuit [23[. Note that a state with a — will not be 
considered a valid input. 

In the following, we show that the output of a post- 
selected quantum circuit can be expressed efficiently as 
a PEPS on a 2D square lattice with both D = d = 2. 
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We start by briefly recalling the concept of measurement 
based quantum computation [Tolflll]: One starts from the 
2D cluster state (which is a PEPS with D = d = 2 [§]) 
and implements the quantum circuit by a sequence of 
projective measurements on the individual spins. Fi- 
nally, the output is found in the unmeasured qubits, up 
to Pauli corrections which depend on the previous mea- 
surement outcomes. In order to express the output of 
a postselected circuit as a PEPS, we therefore start by 
implementing its unitary part in the measurement based 
model. We do this by projecting each qubit on the out- 
come | a) which does not give a Pauli correction, by re- 
placing the original cluster projector Pc with \a)(a\Pc- 
This leaves us with a set of qubits holding the output of 
the circuit, and by projecting the first qubit on |0), we 
obtain the output of the postselected quantum circuit. 
The transformation between the representations can be 
carried out efficiently, and the resulting PEPS has a size 
polynomial in the length of the circuit. 

Conversely, any PEPS can be efficiently created by a 
postselected quantum computer. This holds for PEPS 
on an arbitary graph with degree (the maximum number 
of edges adjacent to a vertex) at most logarithmic in the 
system size, which ensures that the P's are polynomial- 
size matrices. The key point is that any linear map P can 
be implemented deterministically using postsclcction. To 
this end, append rows or colums of zeros to make P a 
square matrix P. By appropriate normalization, we can 
assume w.l.o.g. that P^P < 1. Hence, there exists a uni- 
tary U on the original system and one ancilla such that 
(0|anc^|0)anc = P- This is, by adding an ancilla |0) anc , 
performing U and postselecting the ancilla we can imple- 
ment P. In order to generate a PEPS using postsclcction, 
we thus have to encode each of the virtual spins in [log D~\ 
qubits, create the maximally entangled pairs, and imple- 
ment the C/'s corresponding to the maps P, which can be 
all done efficiently. We are thus left with N ancillas, all 
of which we have to postselect on |0). This, however, can 
be done with a single postselection by computing the OR 
of all ancillas into a new ancilla and postselecting it on 
|0). 

In summary, on the one side we have that any postse- 
lected quantum circuit can be translated efficiently into a 
2D PEPS with D = d = 2, while conversely there is also 
an efficient transform from any PEPS to a postselected 
quantum circuit. In turn, this shows that all the features 
and the full complexity of PEPS can already be found 
in the simplest case of two-dimensional PEPS, making 
them an even more interesting subject for investigations. 



III. THE POWER OF CREATING PEPS 

Let us first briefly introduce the complexity classes #P 
and PP [l2j|. Consider an efficiently computable boolean 
function / : {0,1}^ -> {0,1}, and let s = s(/) := 
\{x : f(x) = 1}| be the number of satisfying assignments. 
Then, finding s defines the counting class #P, while de- 



termining whether s > 2 n ~ 1 (i.e., finding the first bit of 
s) defines the decision class PP. This class contains NP 
and BQP as well as QMA, the quantum version of NP. 

First, we investigate the computational power of cre- 
ating PEPS. More precisely, we consider the scenario of 
Fig. [T] We want to know which decision problems we 
can solve with one use of a PEPS oracle, i.e., a black 
box which creates the quantum state from its classical 
PEPS description, together with efficient classical pre- 
processing and quantum post-processing. 

We now use the PEPS-postselection duality to show 
that the power of creating PEPS equals PP. It has 
been shown that Post BQP — the class of decision prob- 
lems which can be solved by a postselected quantum 
computer — equals PP, PostBQP = PP §. This readily 
implies that a PEPS oracle allows us to solve PP prob- 
lems instantaneously by preparing the output of the post- 
selected circuit as a PEPS and just measuring one output 
qubit in the computational basis. On the other hand, this 
is the best we can do with a single use of the PEPS or- 
acle, since every PEPS can be generated efficiently by 
a postselected quantum computer. BQP postprocessing 
instead of a simple one-qubit measurement does not in- 
crease the computational power, since it commutes with 
the postselection and can thus be incorporated in the 
PEPS. 

The fact that creating PEPS allows to solve PP- 
complete problems strongly suggests the existence of 
PEPS which cannot be created efficiently by a quan- 
tum computer. Note however that the states which ap- 
pear in the PP-hardness proof above are not of this type: 
once the corresponding counting problem is solved, they 
can be easily constructed. While it appears very un- 
likely that all PEPS can be constructed efficiently from 
some normal form (it would imply QMA = QCMA and 
BQP/qpoly = BQP/poly an example of such a state 
is still missing. 



IV. THE CLASSICAL COMPLEXITY OF PEPS 

Let us now investigate the complexity of classically 
simulating PEPS, and its generalization, the contraction 
of tensor networks. For the case of PEPS, there are at 
least three possible definitions of the problem: compute 
the normalization of the PEPS (norm), compute the un- 
normalized expectation value of some observable (uev), 



encoding _ 



BQP post- 
processing 
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no 



FIG. 1: The power of creating PEPS: The original deci- 
sion problem is transformed into a PEPS description by a 
polynomial-time algorithm. The black box creates the corre- 
sponding quantum state, and an efficient quantum postpro- 
cessing returns the solution. Which kind of problems can we 
solve this way? 
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and compute the normalized expectation value (nev). 
Since they can be transformed easily into each other [2J] , 
we will use whichever is most appropriate. 

We first show that contracting PEPS is #P-hard, i.e., 
that for any (polynomial) boolean function /, s(/) can be 
found by simulating a PEPS. Therefore, we take a quan- 
tum circuit which creates ^ x \x) A\f{x))B and encode it 
in a PEPS. Then, the normalized expectation value of <j z 
of B allows to compute s(f). 

To show that the simulation of PEPS is inside #P, 
we have to show that the normalization of the PEPS, 
or equivalently the success probability for the postse- 
lection, can be computed by counting the satisfying as- 
sigments of some boolean function. This can be done 
by adapting well-established quantum complexity tech- 
niques (see Q and references therein): First, approx- 
imate the postselected circuit using only Toffoli and 
Hadamard gates [3, EH • The probability p x for a state 
\x) before postselection is obtained as a kind of path in- 
tegral jl6j |. by summing the amplitudes for all possible 
"computational paths" C = (Ci, ■ ■ ■ jCt-i), where |£ t ) is 
the state step t and T the length of the circuit : 

with a x ,i a product over transition amplitudes Af t _>f t+1 
along the path (. The normalization of the PEPS is 
obtained as the sum over all states where the postselec- 
tion succeeds, J2 x P(o,x)- This can be rewritten as the 
sum over an efficiently computable function /(5, £') = 
a (o.s)x a *o x) (' wn i cn takes values in {0, ±1}, as the cir- 
cuit consisted only of Toffoli and Hadamard gates. Now 
this sum can be computed by counting the satisfying 
assignments of the function /booi(£,z) := (/(£) > z), 
z e {0,1}, which shows that the simulation of PEPS 
is in #P. Together, we find that the classical simulation 
of PEPS is #P-complcte under weakly parsimonious re- 
ductions (see [23]). 

It is natural to ask whether this also shows that con- 
tracting general tensor networks is in #P. For a tensor 
network T, let us denote its contraction by C(T) G C. 
Since the contraction of PEPS is a special case, it is 
clear that the problem is #P-hard. To place it within 
#P, observe first that \C(T)\ 2 = C(T®T*) can be found 
by attaching a physical system of dimension one to each 
site and computing the normalization of the resulting 
PEPS. To determine the phase of C(T), observe that 
C(T © T') = C{T) + C(T'). Thus, by setting T = T* , 
we get |Rc(C(T))|, while the sign can be determined by 
adding another T" = c > 0. This proves that contracting 
tensor networks is #P-complete. 

The obtained hardness results arc stable under ap- 
proximations. To see why, note that any counting prob- 
lem can be reduced to any of our three primitives with 
only linear postprocessing, and thus approximating these 
primitives is as hard as approximating counting prob- 
lems can be. For NEV, this again works by preparing 
Y] \x) A\f{x)) b and computing the expectation value of 



B. For norm and thus UEV, note that the output of 
any normal quantum circuit and thus ^ \x)ji\f(x)) b has 
a known norm when written as a PEPS, since the suc- 
cess probability of each cluster projector is known, and 
the probability of the two measurement outcomes in the 
cluster is unbiased [ll|. Thus, the probability for \1)b 
can be readily determined from the norm of the PEPS 
where we postselected on \1)b- 

V. PEPS AND GROUND STATES 

The interest in MPS and PEPS stems mainly from 
the fact that those states perform extremely well in ap- 
proximating ground states. In the following, we use the 
PEPS-postselection duality, and a relation between post- 
selection and cooling, to shed new light on the connection 
between PEPS and ground states. In particular, we show 
that the unique ground state of a gapped Hamiltonian on 
a /^-dimensional lattice can be approximated efficiently 
by the border of a PEPS with D + 1 dimensions. 

Consider a Hamiltonian on N spins, H = ^;Hi, 
where each Hi acts on a finite number of spins, with 
a unique ground state and a polynomial energy gap 
A > l/poly(iV). Starting from a random state \x), 
the ground state can be efficiently approximated via 
\ipo) ~ exp[— (3H] |x). The imaginary time evolution can 
in turn be approximated using the Trotter decomposi- 
tion, which only requires operations exp[— j3/NHi] act- 
ing on finitely many spins. Since those operations are 
linear, they can be implemented using postselection, and 
we see that postselection can be used to cool into the 
ground state. By embedding the postselected cooling 
procedure in a PEPS, the ground state of any gapped 
iV-particle Hamiltonian can be approximated up to e by 
the boundary of a PEPS, where the extra dimension has 
depth M ~ poly (TV, 1/e) (25[. In case the Hi are local, 
the PEPS can be simplified considerably since any local 
linear operation can be implemented directly on the level 
of the PEPS without the need for ancilla qubits. 



VI. THE POWER OF CREATING GROUND 
STATES 

As we have seen, PEPS can encapsulate problems as 
hard as PP. However, these PEPS are quite artificial, 
while in practice one is often interested in PEPS in con- 
nection with ground states. Therefore, let us have a 
look at the computational power of a ground state oracle, 
i.e., a black box which creates the ground state from the 
Hamiltonian. 

First, let us introduce the complexity class QMA fl7j ] - 
Colloquially, QMA is the quantum version of NP, i.e., 
it contains all decision problems where for the "yes" 
instance, there exists an efficiently checkable quantum 
proof, while there is no proof for any "no" instance. In a 
seminal work, Kitaev [3 has shown that the problem 
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of determining ground state energies of local Hamiltoni- 
ans up to polynomial accuracy is QMA-complete. More 
precisely, in local hamiltonian one is given an N- 
qubit local Hamiltonian H = ^2 Hi with the promise 
that the ground state energy Eq < a or Eq > b, b — a > 
l/poly(7V), and the task is to decide whether Eq < a. 
Clearly, the ground state of H serves as a proof for a 
"yes" instance. In successive works, the class of Hamil- 
tonians has been restricted down to two-particle neareast 
neighbor Hamiltonians on a 2D lattice of qubits [l9[ . 

Let us briefly reconsider our cooling protocol in the 
light of QMA. It is easy to see that the QMA proof need 
not necessarily be the ground state, as long as it is close 
enough in energy (depending on the verifier). Since our 
cooling protocol suppresses higher energy levels expo- 
nentially, the correspondence between postselection and 
cooling shows that a postselected quantum computer can 
be used to create proofs for QMA problems, or differ- 
ently speaking, that any QMA proof can be efficiently 
expressed as a PEPS. 

In the following, we give some observations which indi- 
cate that creating ground states of gapped Hamiltonians 
is easier than creating PEPS. First, note that a ground 
state oracle for arbitrary Hamiltonians is still as power- 
ful as PP. To see why, take a PP problem and encapsu- 
late it in a PEPS. By perturbing the P's randomly by a 
small amount, one obtains a PEPS which is the unique 
ground state of a local Hamiltonian, which can be de- 
rived efficiently from the P's [1, [2(|. This shows that 
an unrestricted ground state oracle enables us to solve 
PP problems. However, the gap A of the above Hamil- 
tonian will be exponentially small: if not, one could add 
a small penalty, say A/100, on the "answer" qubit, and 
use that the original Hamiltonian has ground state en- 
ergy Eq = 0: Then, determining the value of that qubit 
could be solved in QMA, thus proving QMA = PP which 
is considered unlikely 21 1 . 

Since ground states of general Hamiltonians are not 
easier to create than PEPS, let us now assume an oracle 
which only works for local Hamiltonians with a unique 
ground state, known ground state energy, and a poly- 
nomial spectral gap to the first excited state. (Alterna- 
tively, one could consider "proof oracles" for the LOCAL 
hamiltonian problem.) It is easy to see that this re- 
stricted oracle, even with BQP postprocessing, is at most 



as powerful as QMA. The proof is the ground state, and 
the verifier is constructed as follows. Let V\ be the ver- 
ifier for the ground state, it accepts the ground state 
with pqs , and any excited state with probability at most 
Pes = Pes — A, A = l/poly(A r ). Further, let V% be 
the postprocessing circuit which has a polynomial sepa- 
ration between the "yes" and the "no" answer if applied 
to the ground state, p ycs = 1/2 + 5, S = l/poly(iV). 
Take Q — ^ffip , and construct the complete verifier as 
follows: with probability Q, run V\, and with (1 — Q), 
run V2. One can readily check that this gives a polyno- 
mial separation between the cases where the proof is the 
ground state and the postprocessing return "yes" , and 
the cases where either the proof is not the ground state or 
the postprocessing returns "no" . The same strategy can 
be used to show that a PEPS oracle cannot be tested on 
all inputs unless QMA = PP: Otherwise, one could take a 
PP-hard PEPS and construct a verifier which either runs 
the testing routine or reads out the PP solution. 

These observations show that imposing a constraint on 
the spectral gap of a Hamiltonian has direct implications 
on its computational complexity, and we think that the 
complexity properties of gapped Hamiltonians are worth 
being considered. On the one side, in the above scenario 
it is not clear whether all QMA problems can be solved 
using this oracle, on the other side, it is not clear how 
important the knowledge of the ground state energy is — 
note that we however also had this knowledge in the PP- 
hard case. It is also an interesting question whether the 
problem LOCAL hamiltonian remains QMA-complete 
when restricting to polynomially gapped Hamiltonians. 
If not, gapped LOCAL HAMILTONIAN should be a natu- 
ral candidate for a physically motivated class of problems 
weaker than QMA. 
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